Conformal Theories with IR cutoff 
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We give a new perspective on the dynamics of conformal theories realized in the SU(A'") gauge 
theory, when the number of flavors Nf is within the conformal window. Motivated by the RG 
argument on conformal theories with a finite IR cutoff Air, we conjecture that the propagator of a 
meson Gnit) on a lattice behaves at large t as a power-law corrected Yukawa-type decaying form 
GH{t) = ch exp (— mi/t)/t"" instead of the exponentially decaying form CHexp(— mnt), in the 
small quark mass region where mn < cAir: mn is the mass of the ground state hadron in the 
channel H and c is a constant of order 1. The transition between the "conformal region" and the 
"confining region" is a first order transition. Our numerical results verify the predictions for the 
Nf = 7 case and the Nf = 16 case in the SU(3) gauge theory with the fundamental representation. 



Conformal field theories are ubiquitous in nature and 
play important roles not only in particle physics be- 
yond the standard model but also in condensed matter 
physics. Nonperturbative understanding of their dynam- 
ics in four-dimensional space-time is ardently desired. In 
this article, we give a new perspective on the dynamics 
of conformal theories realized in the SU(A^) gauge theory 
when the number of flavors Nf is within the conformal 
window [l[ by studying meson propagators with a finite 
infrared (IR) cutoff. Some preliminary results have been 
presented in [2|. 

Our general argument that follows can be applied to 
any gauge theories with arbitrary representations as long 
as they are in the conformal window, but to be specific, 
we focus on SU(3) gauge theories with Nf fundamen- 
tal fermions ("quarks"). We define the conformal field 
theory in a constructive way We employ the Wilson 
quark action and the standard one-plaquette gauge ac- 
tion on the Euclidean lattice of the size N^ = Ny = N^ = 
N and Nt = rN with aspect ratio r. We impose periodic 
boundary conditions except for an anti-periodic bound- 
ary condition in the time direction for fermion fields. We 
eventually take the continuum limit by sending the lattice 
space a — )> with N oo keeping L = N a fixed. When 
L = finite, the continuum limit defines a theory with an 
IR cutoff. The theory is defined by two parameters; the 
bare coupling constant go and the bare degenerate quark 
mass niQ at ultraviolet (UV) cutoff. We also use, instead 
of go and mo, /3 = G/^g and K = l/2(moa-|-4). For a later 
purpose, we define the quark mass ruq through Ward- 
Takahashi identities with renormalization constants be- 
ing suppressed [2[. 

Let us quickly remind ourselves of the renormalization 
group (RG) flow when Nf is in the conformal window. 
One important fact is finite size lattices in computer sim- 
ulation always introduce an IR cutoff Air ^ l/{Na). If 
the IR cutoff were zero, when quarks have tiny masses, 
the RG trajectory would stay close to the critical line, 
approaching the IR fixed point and finally would pass 
away from the IR fixed point to infinity. Therefore the 
IR behavior is governed by the "confining region" . Only 



on the massless quark line the scale invariance is realized 
at the IR fixed point. See the left panel of FIG. 1. 

When the cutoff Air is finite, the RG fiow from UV 
to IR does stop evolving at the scale Air. When the 
typical mass scale (e.g. that of a meson) uih is smaller 
than Aji^, it is in the "conformal region". On the other 
hand, when mn is larger than Air, the flow passes away 
from the IR fixed point to infinity with relevant variables 
integrated out, thus being in the "confining region". See 
the right panel of FIG. 1. 

This scenario implies that when physical quantities at 
IR (e.g. hadron masses) arc mapped into a diagram in 
terms of physical parameters at UV (e.g. the bare cou- 
pfing constant and the bare quark mass), there will be 
gaps in the physical quantities along the boundary be- 
tween the two phases. There the phase transition will be 
a first order transition. 

To make the above RG argument more concrete, we 
will study the propagator of the local meson operator 



(1) 



When the theory is in the "confining region", it decays 
exponentially at large t as 



Gnit) = CH cxp(-tohO 



(2) 



due to the physical one-particle pole, where mn is the 
mass of the ground state hadron in the channel H. 

In contrast, we claim that in the "conformal region" 
defined by 



TTT-H < cA 



IR, 



(3) 



where c is a constant of order 1 which we will determine, 
the propagator G{t) behaves at large t as 



GH{t) 



exp(— m^ft) 
^« ^ ' 



(4) 



which is a power-law corrected Yukawa-type decaying 
form instead of the exponential decaying form (cq.(P|)) 
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observed in the "confining region" . We further ciaim 
that the boundary between the "conformai region" and 
the "confining region" is a first order transition. 

We note that the behavior eq.(|3]) is proposed based 
on the AdS/CFT correspondence with a softwall cutoff 
in the iiterature The meson propagator in the mo- 
mentum space has a cut instead of a poie: Gh{p) = 
The propagator in the position space 
(after space integration) tak;cs the form eq.Q in the limit 
trfiH > 1. 

We distinguish ttih in eq.Q from the pole mass mn in 
eq.®. Eq.Q thus means the lower limit of itih is c Air . 

In the continuum limit with L = oo (i.e. Air = 0), the 
propagator on the massless quark line takes the form 

GH{t)^c^, (5) 

consistent with itih = limit of eq. (jj]). If we take the 
coupling constant go = g* at the UV cutoff, an takes a 
constant value, and the RG equation demands 

«if = 3-27*, (6) 

for the pseudo-scalar (PS) channel with 7* being the 
anomalous mass dimension ^ at g — g*. The theory 
is scale invariant (and shown to be conformai invari- 
ant within perturbation theory See also e.g. Q 
and references therein from AdS/CFT approach). When 
< 5o < ff*: Q^ff depends slowly on t as a solution of 
the RG equation. In the IR limit f — > 00. we must retain 
aH{t) ^ 3 - . 

The continuum limit with L ~ finite defines a con- 
tinuum theory on X which corresponds to a com- 
pact three-torus at finite temperature. The IR cutoff 
Air ^ 1/i is finite. The propagator Gnit) behaves at 
large t as a power-law corrected Yukawa type decaying 
form eq.([3]). The exponent an in t 00 with tfriH <C 1 
takes the universal formula 3 — 27* while with t rfiH ^ 1 
it takes a value depending on the dynamics, from which 
we can estimate the mass anomalous dimensions (see [2[ ) . 

Now let us discuss the results of our numerical calcula- 
tions. We have performed simulations in the Nf = 7 and 
Nf = 16 cases, which we conjecture Q, are the bound- 
aries of the conformai window. The algorithms we em- 
ploy are the blocked HMC algorithm for Nf = 2N 



and the RHMC algorithm g for Nf = 1 in the case 
Nf ^2N + 1. 

We specify the coupling constant /? = 6.0 for Nf = 7 
and (3 = 11.5 for Nf = 16, taking account of the fact that 
the IR fixed point for Nf = 16 is /3 = 11.48 in two- loop 
approximation. We use the lattices of size 16^ x 64 and 
24 X 96, and vary the hopping parameter K so that the 
quark mass takes the value from 0.40 to 0.0: We simulate 
with 15 hopping parameters on the 16'^ x 64 lattices, and 
5 hopping parameters on the 24^ x 96 lattices. 

We choose the run-parameters in such a way that the 
acceptance of the global metropolis test is about 70%. 
The statistics are 1,000 MD trajectories for thermaliza- 
tion and 1,000 MD trajectories or 500 MD trajectories 
for the measurement. We estimate the errors by the jack- 
knife method with a bin size corresponding to 100 HMC 
trajectories. 

Let us first discuss the results for the Nf — 7 case on 
the 16'^ X 64 lattices. We define the effective mass mH{t) 

bv c°sh(mH(t)(t-Aft/2)) _ Gait) j^^p r, „u„^„ f 

dependence of the effective mass for the PS channel with 
three types of sources. On the left panel, we see the clear 
plateau of the effective mass at t = 22 '--^ 31 when quark 
mass is relatively large; ruq = 0.25 (A' — 0.1400). On 
the other hand, on the right panel, we see the effective 
mass is slowly decreasing without no plateau up to t = 31 
when the quark mass is small; ruq = 0.045 [K = 0.1459), 
suggesting the power-law correction. 

The effective masses for all cases with rriq < 0.172 
(A' > 0.1415) exhibit a similar behavior to the K = 
0.1459 case. We show the power-law corrected fit for the 
local- local data in the K — 0.1459 case with the fitting 
range t ^ [15 : 31] on the left panel in FIG. 3. The fit 
with an = 0.8(1) reproduces the date very well. 

In the case of Nf = 16, the effective masses on the 
16^^ X 64 lattices exhibit the power-law corrected Yukawa- 
type decaying form for all cases when < 0.0978 (AT > 
0.130). We show the power-law corrected fit for the local- 
local data in the K = 0.1315 case with the fitting range 
t = [15 : 31] on the right panel in FIG. 3. The fit with 
aH = 1.15(2) reproduces the data well. 

We also report that the data on 24^ x 96 lattices for 
both cases of Nf = 7 with K > 0.1459 and Nf = 16 
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FIG. 2. The eflFective mass for Nf = 7 (left: for K = 0.1400; exponential decay) and (right: for K — 0.1459; power-law corrected 
Yukawa decay): With three types of sources; the local-sink local-source (black squares), local-sink doubly-smeared-source (red 
circles) and local-sink doubly- wall-source (green triangles). 
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FIG. 3. The effective mass plots for local-sink local-source case and fits by power-law corrected Yukawa type decay (left: for 
Nf = 7 and K = 0.1459) and (right: for Nf = 16 and K = 0.1315). 



with K > 0.130 clearly exhibit the power-law corrected 
Yukawa type decaying form. We obtain an = 0.54(17) 
for K = 0.1459 in iV/ = 7 and an = 1.57(28) for 
K = 0.130 in Nf = 16, respectively. Thus we have con- 
firmed that the propagators behaves as the power-law 
corrected Yukawa type decaying form, when the quark 
mass is small enough. 

In the remaining part of the article, we try to iden- 
tify the boundary between the "conformal region" and 
the "confining region" more precisely by investigating the 
transition region carefully with a small step of the value 
of K for both of the Nf = 7 and Nf = 16 cases. 

We present the results for and mps (or rhps) in the 
Nf = 7 case in FIG. 4. We first note that the quark mass 
•niq denoted by black points and line on the left panel is 
excellently proportional to 1/ K in the whole region from 
0.006 to 0.555. 

For the propagators of the PS meson, we observe the 
clear transition from the exponentially decaying form to 
the power-law corrected Yukawa-type decaying form at 
K = 0.1412 - 0.1413. For K < 0.141, we present mps 



obtained from the fit to the exponentially decaying form 
with the fitting region [28:31] by green points and line. 
For K > 0.1415, we present rhps obtained from the fit 
to the power-law corrected Yukawa type decaying form 
with fitting range [15:31] by blue points and line. 

The transition region is enlarged on the right panel in 
FIG. 4. The values of mps and rhps are different in the 
limit K = 0.1412 from smaller K and from larger K. It 
should be noted that even if we had ignored the power- 
law corrections in Yukawa type decay and had estimated 
the nips (plotted in FIG. 4 by black points) by assuming 
exponential decay, there would have been a gap at the 
transition point. 

It is worthwhile to stress that at K ~ 0.1412 and 
0.1413 we observe two states, depending on the initial 
sate. The existence of two states persists at least with 
500 MD trajectories. We also note that the exponent 
an does not vanish toward the transition point; it rather 
increases up to aa = 1-1 ^ 1-3. Thus the transition is 
a first order transition at K = 0.1412 ^ 0.1413. These 
facts imply that there indeed exists a "conformal region" 
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FIG. 4. rriq and mps (or rhps) vs. 1/A' for the range 0.130 < K < 0.1472. The transition region is enlarged on the right 
panel. See the text for the symbols. 



which is separated by a gap of physical quantities from 
the "confining region" at finite quark mass. 

The TTiq dependence of rhps is rather comphcated. 
Once rhps becomes small at the transition point, it 
slightly increases as nig decreases, and then decreases 
at K < 0.143. Apparently the small "rhps region suffers 
from finite size effects. To verify the scaling relation [l^ 
for fnps in terms of niq we have to control finite size 
effects. 

Taking mps = 0.601 - 0.625 with ruq = 0.177 - 0.183 

{K = 0.1413 ^ 0.1412), as the critical mass in eq.(3), we 

estimate c = 2.16 ^ 2.25, with our working definition of 
Air, = 2 7r(Af3 

In the Nf = 16 case we have observed similar results. 
The niq dependence of rhps is similarly complicated. The 
transition occurs at mps = 0.513 ^ 0.539 with niq = 
0.237 - 0.244 {K = 0.1255 ~ 0.125) from which we 
estimate c = 1.85 ~ 1.94. 

In the cases of 24"^ x 96 lattices for both Nf = 7 and 
Nf = 16 cases, we have also observed the transition : In 
Nf = 7, we find exponential decay for niq > 0.061 {K < 
0.1452) and power-law corrected Yukawa type decay for 
niq < 0.045 {K > 0.1459). In Nf = 16, we find power- 
law corrected Yukawa type decay for niq < 0.0978 {K > 



0.130). The critical mass decreases compared with the 
16'^ X 64 lattices as expected since the IR cutoff decreases. 

We conclude the paper with two main goals achieved 
in our study. We verified (1): The meson propagator 
shows a transition from an exponentially decaying form 
to a power-law corrected Yukawa-type decaying form at 
the critical hadron mass whose value is specified by the 
IR cutoff. (2): The transition is first order. 

In an accompanying publication [T]| . we will discuss 
the physical interpretation of the power-law corrections 
in terms of the dynamics of the underlying conformal 
field theory by taking the continuum limit. We also hope 
to report the physical applications of our findings in the 
near future. 
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